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Abstract 

The goal of this article is that of understanding how the oscillation 
and concentration effects developed by a sequence of functions in M"^ 
are modified by the action of Sampling and Reconstruction operators 
on regular grids. Our analysis is performed in terms of Wigner and 
defect measures, which provide a quantitative description of the high 
frequency behavior of bounded sequences in (R'') . We actually 
present explicit formulas that make possible to compute such measures 
for sampled/reconstructed sequences. As a consequence, we are able 
to characterize sampling and reconstruction operators that preserve 
or filter the high-frequency behavior of specific classes of sequences. 
The proofs of our results rely on the construction and manipulation of 
Wigner measures associated to sequences of discrete functions. 
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1 Introduction 

1.1 Statement of the problem: oscillation and concentration 
under the effect of sampling and reconstruction 

A central problem in Numerical Analysis and Signal Theory is that of recon- 
structing a function u (x) defined in from a discrete set of measurements 
taken on an uniform grid of step size h. This discrete values are typically 
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obtained by applying to the function u a sampling operator 5^ of the 
following type: 

for some sampling function (p. One then tries to recover u by means of a 
reconstruction (or interpolation) operator through the formula: 

T^^5> (x) := ^> (^) - ^) > (1) 

where ijj is some fixed reconstruction function. This process usually only 
provides an approximation of the original function u, with an error that 
vanishes as h tends to zero. Such reconstruction schemes have been the 
object of intensive study both from the point of view of Approximation 
Theory and Numerical Analysis. 

Here we shall be concerned with the high-frequency approximation 
properties of those operators, that is, we shall study how a reconstruction 
scheme such as is able to capture (or filter) oscillation and concentra- 
tion-like phenomena on the functions it is intended to approximate. More 
generally, me shall be interested in clarifying how the high frequency be- 
havior of a sequence of reconstructed functions depends on the profiles ip 
and the sampling rate h chosen. 

Before giving a more precise statement of our objectives, let us first 
illustrate the above discussion with two specific examples: consider (x) := 
k'^/'^p {k {x — xo)) and (x) := p (x) e^^^'^ with p G (M'^) ; the sequence 
(/fc) concentrates around the point xq as A; ^ oo, whereas {gk) oscillates 
in the direction . The results we shall present in this paper are aimed 
to understand to what extent the sequences {r^p^ S^'^ f]^ and (t^'^ Sj^^ gk^ 

reproduce the same behavior as (fk) and (gk) (i-e., if concentration and 
oscillation persist), for a given sequence (hk) of positive reals that tends to 
zero (the sampling steps) and some choice of 99 and ip. 

Perhaps, the simplest convenient setting to formulate our results is pro- 
vided by the notion of defect measure, an object that gives a quantitative 
description of what we shall understand by concentration and oscillation 
effects and whose definition we next recall. Let (n^) be a weakly converg- 
ing sequence in the space (M*^) ; denote by u its weak limit and remark 
that the densities \uk — are uniformly bounded in (M*^). Helly's com- 
pactness Theorem then ensures that some subsequence ( \uk„ — u\'^] weakly 
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converges in the set of positive Radon measures; or, in other words, that 
there exists a positive measure v on M"^ such that 

/ (/> (x) (x) — ti (x)l^ ^ / 4>(x)dv(x) as n ^ cxo, 

for every (j) ^ Cc (W^^ . When the above convergence takes place without ex- 
tracting a subsequence we say that v is the defect measure of the sequence 

Immediately from this definition one deduces the following general prin- 
ciple: if V is the defect measure of a sequence {uk) and a; C M*^ is a bounded 
Borel set, then there is an equivalence between v [uj) = and the fact that 
tifclo; converges strongly to in {uj). Thus, the support of v is precisely 
the set where strong convergence fails, that is, the set where oscillations and 
concentrations take place. 

But defect measures are also able to detect concentration and oscillatory 
phenomena and give quantitative information about them. Consider the 
sequences (/fc), {gk) previously defined; they both weakly converge to zero 
in (R^) and it is easy to check that their respective defect measures 
are ||/o||^2(']Rd^ <^a'o ^'^'^ \p{x)\^ dx. Notice that, in the first case, the defect 
measure actually captures the concentration of the sequence around the 
point X = xq. In the complementary of that point, where the sequence 
converges strongly to zero, the measure vanishes. In the second example, 
the defect measure is uniformly distributed on M*^, this being consistent with 
the fact that strong convergence does not take place in any subset of W^. 

Let us point out that the analysis of concentration and oscillation effects 
developed by a sequence of functions is a central issue in many problems of 
the Calculus of Variations and Partial Differential Equations. A number of 
applications of defect measures may be found in the analysis of variational 
problems with loss of compactness performed by P.-L. Lions in 112j .^ 

Consider a sequence (wfc), weakly converging to zero in (R'^) ; sample 
it using a profile if and form the reconstructed sequence 

^From now on, we shall use the term measure as an abbreviation of the longer Radon 
measure. Recall that the space of Radon measures M (K'') is identified, by Riesz's Theo- 
rem, with the space of continuous linear functionals on Cc (R'') . 

^We also refer to L.C. Evans' notes 0] for an exposition of some additional applications 
as well as a discussion of other measure-theoretical objects (such as, for example. Young 
measures) designed to study the failure of strong convergence. 
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for some given ip and some sequence (/ifc) of positive reals tending to zero. 
The functions Vk are bounded in (M*^) and tend weakly to zero provided 
ip and ip satisfy suitable hypotheses (see Lemma |3.1l in Section |21 below) . 
Suppose furthermore that the densities \vk\ weakly converge to the defect 
measure i^^^^. 

One of the main issues addressed in this article is that of understanding 
the relations existing between the defect measure z^,^,,/,, the profiles ip, ij) and 
the sequences (n^), (/ifc). Among these, we point out: 

A. Is there a formula, valid for any sequence (n^), relating v^^^ to the 
defect measure v only in terms of the profiles 99 and V'? 

B. Given (n^), characterize the profiles 99 and such that z^^^^ = 0. 
This is the problem of filtering since, as we have discussed before, z^^^^ = 

is equivalent to the strong convergence to zero of the sequence ^T^* 5*^* 

C. Similarly, characterize the profiles ip and ip such that i^<^^^ = v for a 
given {uk). 

D. Finally, characterize the profiles that give v^p^^ = v for every (u^). 

We shall prove that the answer to question A is negative. This is due 
to the fact that the measure i/^^^ is sensitive to the characteristic directions 
of oscillation of the sequence (u^), whereas v is unable to distinguish them. 
As we have seen above, the defect measure of the oscillating sequence (5^) 
equals |/9 independently of the vector ^\ that is not the case for v^p^^. 

Indeed, under additional assumptions on Lp and ip we prove (see Theorem 
Oand Corollary OJ: 



V 



(^)= E |V^(^° + 2vrA;) \p {e)Up{xt dx. 



k& 



Thus the measure i^^^^ is .^'^-dependent and cannot be expressed solely in 

terms of z^, 99 and ip. Note that z^^^^ is identically zero as soon as any of 

^ 2 

X^fcezd V' {^^ + 27rA;) or ip is null. Analogously, the profiles that give 
^tfi,ip = ^ are precisely those which satisfy 



Therefore, in order to understand how u^^^ is built, we must have at our 
disposal an object that is able to distinguish between oscillatory phenomena 
at different directions. 
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1.2 Wigner measures 

This refinement is provided by the theory of Wigner measures.^ Given 
a bounded sequence in one associates to it a measure on 

X M'^ which describes the concentration and oscillation effects (these are 
the respective roles of the variables x and ^) occurring at some characteristic 
length-scale. This measure takes into account the characteristic speeds as 
well as the directions of propagation of oscillations. One way of defining 
them consists in replacing the density involved in the definition of 

the defect measure by the phase space (microfocal) density: 

M (x, := i^e) e-«/% (2) 

where u is the Fourier transform of u and e is a positive constant. The 
{2tt)~'^ factor in the definition of m'^ [u] is placed to have: 

m^u]{x,OdC = \u{xf, [ m^u] (x, g) dx = '^/^^^f . (3) 

Thus, the function [u] may be looked at as joint physical space-Fourier 
space "density" , in spite of the fact that [u] is not positive in general. 
However, limits of these quantities are positive measures: 

Theorem 1.1 Let (u^) be a bounded sequence in L^ (M'*) and let (e^) be a 
sequence of positive numbers tending to zero. Then it is possible to extract 
a subsequence (ufc„) such that, for every test function a £ S {R'^ x W^) , 

lim / a{x,(,)m'"'^ [ukJix,C)dxd£, = a {x,(,) dfi {x,^,) , (4) 



where ji is a finite positive measure on x M*^. 

A measure // € M.+ {W^ x M*^) is called the Wigner measure of the 
sequence (u^) at scale (e^) whenever the limit Q) holds without extracting 
a subsequence. Different proofs of Theorem 11.11 mav be found in [51 IT^ 17]. 
Let us point out that other quadratic densities may used to define Wigner 

^This object is present in the work of E.P. Wigner on semiclassical quantum mechanics 
im. Recently, Wi gner measures have gained interest since the works of P. Gerard 0, P.-L. 
Lions & Th. Paul [HI, P. Markowich, N. Mauser & F. Poupaud |T1| among others. Related 
objects are the Microlocal defect measures or H-measures, introduced independently 
by P. Gerard |Sj and L. Tartar [TH]. 
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measures. For instance, in // is obtained by replacing rrf [u] in the limit 
Q, by the more familiar Wigner transform: 

dp 



[u] {x,i) :-- 



u X 



u { X + e 



P 



Jp<. 



(2^) 



(5) 



It is also possible to consider Wave-packet (Husimi) transforms. Of 

course, all this methods are equivalent (the same limit is obtained), cf. the 
discussion in [S]. 

The Wigner measure encodes all the information contained in the defect 
measure provided the sequence (uk) oscillates at frequencies of the order of 
e^^. More precisely (see [51 IT^): 

Proposition 1.2 // fi is the Wigner measure at scale (e^) of a sequence 
(itfc) and V is the measure obtained as the weak limit in J^j^ (M'^) of the 

2 

densities \uk\ dx, then the identity 

v{x)= I fj. (x, dC) 
holds provided (u^) is e^- oscillatory: 



limsup / \uj:{0\^d^ 







as R 



oo. 



(6) 



Notice that condition © actually expresses that the energy of the Fourier 
transform of Uk is concentrated in a ball of radius R/sk, which should be 
understood as the requirement that the sequence (li^) does not oscillate at 
length scales finer than e^. 

To illustrate this discussion it may be helpful to look at explicit com- 
putations. The Wigner measure at scale (ek) of the concentrating sequence 
(/fc) defined at the beginning of this section is given by: 



/i(x,^) 



I l|2 

\p\\l2 



Sxo {x) (g) 6o (0 if Skk 



0, 



JL 

(2vr)' 



if Sk 



(7) 



if Ekk oo, 

while for the oscillating sequence {g^) it can be checked to be: 



/i(x,0 



|p (x)l^ dx (g) (0 if efe/c ^ 0, 



|/>(x)| dx0 6^o{S^) iiek = k 
if Ekk 



-1 



(8) 



oo. 
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These examples show the importance of the choice of the scale (efc). When 
this scale is taken to be coarser than the characteristic length-scale of 
oscillation/concentration, it is no longer true that the projection on the first 
component of their Wigner measures coincides with the defect measure. On 
the other hand, in the case e^k (the scale chosen is much smaller than 
the actual oscillation scale) the Wigner measure is not able to capture the 
direction of oscillation. Hence, to obtain a complete description, the scale 
(efc) must be taken of the same order than that of the oscillations. 

Wigner measures turn out to be the correct tool for comparing the high 
frequency behavior of the sequences {uk) and (Tj^'^Sll'' 



1.3 Computation of Wigner and defect measures 

Given a sequence of sampling steps (/ifc), it is clear that the functions 
T^^ S^'' Uk will not develop oscillation and concentration effects of character- 
istic sizes asymptotically smaller that h^. Most commonly, these functions 
will form an /ifc-oscillatory sequence;^ consequently, only Wigner measure at 
scales coarser or of the same order than (/i^) will be considered. 

In order to establish explicit formulas, we shall require additional hy- 
potheses on If, and on the Wigner measures involved. Nevertheless, in 
order to simplify the statement of our results, in this introduction we shall 
impose the following (more restrictive) condition on the admissible profiles: 



It {x)\<C{l + Ixl)"''"^ , for every xGR"^ and some C,e > 0. (9) 

More general results may be found in Section [3 
We prove the following: 

Theorem 1.3 Let ip, tp satisfy Suppose (uk) is a bounded sequence in 
(W^) and that /i is its Wigner measure at scale (hk). Suppose moreover 
that the measures 

\^p{^ + 2^Tn)\^ fi{x,^ + 2TTn) (10) 

are mutually singular for n £ Z^. 

Then the Wigner measure at scale (hk) of the sequence {t^^ S^^'^Uk^ is 
given by: 

— I ^ — ^ 2 

/^<^,v(a^'0= '4'{i)\ 2^ l<?(C + 2vrn)| /i (x, ^ -h 27rn) . 



'However, this may fail for some pathological examples (see paragraph 15. 31 1. 
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From this, one deduces: 



Corollary 1.4 //, moreover, 
v^^^ then: 



'lj*S^''U}^ dx weakly converges to a measure 



{x)= I |v.(^ + 27rA:) ' 1^(01 



This shows, in particular, that a formula relating and the weak 
limit V of \uk\^ dx does not exist unless {uk) is /ifc-oscillatory and is 
of the form v [x) ® It also shows that v = Vip^^ if and only if 



there do not exist profiles V satisfying © such that v equals i^^^^ for 
every /ifc-oscillatory sequence (ufc). 

On the other hand, Theorem ll.3l imDlies that question A above does have 
a positive answer in terms of Wigner measures, at least when restricted to 
the class of sequences which satisfy (jlUp. That condition, roughly speaking, 
imposes a restriction on the size of the region in frequency space where 
an admissible sequence fails to converge strongly to zero. Below, we shall 
compare it with that appearing in Shannon's sampling Theorem. 

The above results will be obtained as corollaries of the more general 
Theorems 17. II and 17.31 Profiles that belong to negative-order Sobolev spaces 
or that fail to satisfy the localization hypothesis Q are allowed. However, 
this will require to impose compatibility conditions on the Wigner measure 
II. 

As an illustration of the range of results that will be obtained in this 
more general setting, we present an asymptotic version of Shannon's 
sampling Theorem.^ It corresponds to taking as sampling profile ip = 5o, 
the Dirac delta at the origin, and as reconstruction function tp := 1q, where 
Q := [— vr, tt)'^. Notice that S^^u (n) = u (hn) is the discretization operator, 
whereas the corresponds to band-limited reconstruction. 

Theorem 1.5 Let (uk) be a bounded sequence in (M'^) and denote by fi 
its Wigner measure at scale {hk). Suppose, in addition, that G (M'^) 
for some s > d/2 and 



YlkeZ'i V + 27rA;) IP (01^ = 1 for /i-almost every ^ G M.'^. Consequently, 




(11) 



^See paragraph 13. 1 1 for a statement of Shannon's original sampling Theorem. 
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Then, the Wigner measure at scale {hk) of (t^J^'' S'^I^Uk^ is 



(12) 



Moreover, if 



spectively, then 



dx and \uk\ dx weakly converge to CLnd v, re- 



{x)= I M (x, di) = v{x) . 



Thus, unlike the operators considered in Theorem ll.31 the composition of 
discretization and band-hmited reconstruction preserves the defect measure 
for a large class of sequences. 

Notice that, by the Sobolev imbedding Theorem, S^^u^ is well-defined. 
Actually, (jTTl i^ ensures that the sequence of discretizations is square-summable 
and, consequently, that {t^'' S^^^Uk^ is bounded in (M*^) and /ifc-oscillatory 
(for a more complete result, we refer to Lemma l3.1() . Condition (|lll ii) ap- 
pears because il) is not continuous; we shall discuss its necessity in paragraph 
14.41 Finally, (llll iii) should be understood as the analog of Shannon's original 
band-limited condition in this context. 

To conclude this short description, let us present how the above results 
may be refined when the sequence {uk) is known to be efc-oscillatory and 
the sampling rate {hk) is taken to satisfy hk/sk 0. As it can be expected, 
much more precision is gained: 

Theorem 1.6 Suppose ip satisfy ^ and {uk) is an £k-osciUatory, bounded 
sequence in (M*^). If fi is its Wigner measure at scale {sk) then the cor- 
responding measure of the sequence ( T^'^S'^'^Ufe j is 



l^ip,4> 



V(0) |i?(0)|V- 



Moreover, if the densities a 
^ip,^) o.iT'd V respectively then 



T^'Sl^'^Uk 



dx and \uk\ dx weakly converge to 



E\ — 2 
V'(27rn) 1^(0)1 u{ 
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This Theorem holds under much more general conditions on (/? and ip 
(see Theorem 17. 6() and gives a positive answer to question A provided we 
consider only e^-oscillatory sequences. 

An immediate consequence of the above result is that zero-mean sam- 
pling profiles (p (i.e. with (p (0) = 0, as a wavelet, for instance) completely fil- 
ter any oscillations that occur at scales much coarser than the sampling rate 
hk- For such a profile, u^^^ = for every Sfc-oscillatory sequence. An anal- 
ogous phenomenon occurs for reconstruction profiles satisfying ip (27rn) = 
for every n G Z*^. 

On the other hand, a sufficient condition to have equality between z^^ ,^ 



and 1/ is that |^ (0)| = ^ (0) =1 and (27rn) 



for n 7^ 0. 



1.4 Strategy of proof: Wigner measures in the discrete set- 
ting 

The proof of the results we have presented above will be achieved by ana- 
lyzing separately the sampling and reconstruction operators Sj^ and T^. In 
order to develop, it is necessary to deal with a concept of Wigner measure 
associated to a sequence of discrete functions. We shall introduce it 
by means of a discrete analogous of the transform [■]. We detail this in 
the following paragraph. 

To a discrete square-summable function U ^ L"^ (/iZ"^), where (/iZ'^) 
stands for the space of the functions U defined on Z'^ with values in C such 
that the norm 

is finite, we associate: 

y.2d ^ /h \ 

[U] {x, := ^rnU e^"^-(V-)«5,„ (x) . (13) 

Here, 6hm is the Dirac mass centered at the point hm and U denotes the 
discrete Fourier transform: 



which, as is well known, is a 27rZ'^-periodic function in Lf^^ (R*^) . The dis- 
crete transform [U] may be related to the continuous [u] by noticing 
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that 

where T^U (x) = h'' U^^dhn (x) . (14) 

This is meaningful, since [u] is weh-defined for any tempered distribution 

In order to simphfy our language we make the following definition: 

Definition 1.7 Let h = (hk) be a scale. We shall call a sequence [U^''^ 
hk'bounded if and only if U^'' G L?' [hkl-''') and ^U^''\\y^ < C for every 
/c G N. 



[U] 



m 



One has the following convergence result (which is not a direct conse- 
quence of Theorem ll.ip : 

Proposition 1.8 Let (hk), (sk) scales such that {hk/£k) bounded and 
let (U^f'^ be an h^-bounded sequence of discrete functions. Then i^M^'^ [C/'^'^]) 
is hounded in S' (W^ x M*^) and given any of its convergent subsequences 
(U^kn\ there exists a positive measure /i such that, 



lim ( M^'^" 

n^oo 



,a) = / a{x,C)dfj,{x,(,) , (15) 



for every a £ S 

This will be proved as a Corollary of the more general Proposition 13.41 
which in turn follows from the analysis of Wigner measures in negative- 
order Sobolev spaces that is performed in Section |HJ As in the continuous 
setting, we say that a measure /u is the Wigner measure at scale (e^) of 
a sequence of discrete functions if the limit H15|) holds for the whole 

sequence. 

Remark 1.9 i) When (hk/sk) is unbounded, it may happen that M''''" 
is not bounded in S' (M*^ x R'^) . 

a) If hk/sk ^ c > then fi is not finite. Indeed, it is periodic (with 
respect to the lattice (27r/c)Z'^J in the ^ variable. 

Hi) However, when hk/sk 0, the Wigner measure /i is finite, as in the 
continuous case. 
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With this tool at our disposal, we are able to compare the Wigner mea- 
sure of a sequence of discrete functions ([7^*=) with that of a reconstructed 

sequence (^T^'° U^'' ^ . Analogously, we may compute the Wigner measures 

of sequences of sampled discrete functions i^S^'^Uk) in terms of those corre- 
sponding to the original sequence (itfc)- These are respectively the contents 
of Theorems 14.61 and 14.21 

1.5 Plan of the article 

Results and assumptions concerning the operators 5^ and are collected 
in Section 131 

In Section 12 the problem of computing Wigner measures for sequences 
of sampled or reconstructed functions is addressed. Formulas for Wigner 
measures at scales of the same order than the sampling/reconstruction step 
(/ifc) are presented in Theorems 14.51 and 14.21 Theorems 11.31 and 11.51 then 
easily follow from those two results. We also point out the relationships 
existing between these Wigner measures and the concept of Wigner series 
introduced in |^ . 

The problem of the computation of defect measures of sequences of the 
from {rlp' U^''^ is considered in Section [SJ the main results are presented in 
Proposition 15.81 and Corollary 15.91 

In Section El we investigate Wigner measures at scales (e^) satisfying 
hk/^k ~^ 0. Explicit formulas are presented in Theorems 16.11 and 16.21 from 
which Theorem II .61 immediatelv follows. 

The composition of sampling and reconstruction is studied in Section [71 
the main results of this article are proven there. 

Finally, Section jSl contains the elements from the Theory of Wigner mea- 
sures on which the proofs of most of the results of this article are based 
on. Propositions 18. ll and 18.31 which extend the Theory of Wigner measures 
to sequences in Sobolev spaces of negative order, are systematically used 
throughout this paper. 

2 Notations and conventions 

We briefly present some notation that will be used throughout this article. 

B (x; R) will denote the open ball with radius R of centered at the 
point X. We shall set 

Q ■= [-7r,7r)'^, 
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and 1a will denote the characteristic function of a set A C W^. 

We write F to denote dc lattice 27rZ'^. A function / defined on is 
F-periodic if / (x + 7) = / (x) for every 7 G F and every x G M"^. 

We adopt the following convention for the Fourier transform: 



n(^) := [ u{x)e-'^-^dx. 



Given a measurable function (/^(C), the Fourier multiplier of symbol 
if is the operator {D^) formally defined by 

^{D,)u{x):= [ ^i^)u{Oe''<-^=^*u{x), 

(p being the inverse Fourier transform of Lp. 

A particularly important Fourier multiplier is the Bessel potential 
(Da;), of symbol 

(0 := fi + ier 



Next, we recall the definition of some function spaces. 

As usual, S (M'^) denotes the space of rapidly decreasing functions 
and S' (R'^) stands for its dual, the space of tempered distributions. 

Given r G R, W (R"^) , the Sobolev space of order r, consists of the 
distributions tt G 5' (R'^) such that {D^Y u e (R''). 

The weighted space (R'^; {xf) is that of the functions u G L\^^ {W^) 
such that 

\ 1/2 

l^tll r2('ra)d./„\r\ := I / \U (xW (x)' CLX \ < OO. 



/ \u{x)\' {xY dx] 
Jm^ J 



The analogous definition is understood for (R*^; {xY)- 

By C°° [W^; (xY) we intend the space of functions u G (R'^) such 
that 

ll^"'"llL°°(]Rrf-(a;)'') ^ °° every multiindex a G N''. 



Co (R ) denotes the spaces of continuous functions on R" vanishing at 
infinity. 

Given an open set C R'^, M.+ (O.) is the set of positive Radon 
measures on Q,, which can be identified through Riesz's Theorem to the 
set of positive functionals on Cc {^) , the space of continuous functions on J7 
with compact support. 

Finally, in order to lighten our writing, 
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we 



shall write S and S' instead of S {RixRi) and S' [Rix 




For a measurable function / : C, we use the notation 

Df := <x eM."^ : / is not continuous at x > . 



An important, perhaps non-standard, definition is that of a scale: 

Definition 2.1 A scale (efe) is a sequence of positive numbers that tends 
to zero as k ^ oo. 

Given two scales {hk) and (sk), the notations hk <C ejt and hk ~ £jfc 
will be used to indicate that limk_,^hk/ek = and limjt^oo /ife/e^ = c > 
respectively. 



3 Sampling and reconstruction 
3.1 Definitions and examples 

The sampling and reconstruction operators wc arc going to consider are next 
described. Given a distribution tp e S' (W^) we set for every n G Z'' and 
h>0, 



The reconstruction (or synthesis) operator T^, acting on discrete func- 
tions U of Z'^ is defined to be 



This expression is well-defined for finitely supported discrete functions. When 
is a continuous function such that ip (0) = 1 and 99 (k) = for /c G Z'^ \ {0} 
then T^^U is actually a function that interpolates the discrete values Un 
on the grid /iZ'^, i.e. T^U (hn) = Un for all n G Z'^. 

Analogously, the sampling (or analysis) operator 5^, a priori only 
acting on functions u e S (W^) , is defined as follows: Sj^u is the discrete 
function given by 





(16) 
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When Lp = (5o, we obtain the usual discretization operator: S^^u{n) = 
h~'^u [hn) for every n G TJ^. 

Indeed, these sampling/reconstruction schemes include several well-known 
such procedures on regular grids. Among many others we may cite: 

• Cardinal B-Splines. The S-spline of order zero is the function 
(x) := lj_-^^2 1/2]'* (^)' function TJ^U is just the piecewise constant 

interpolation of the discrete function U on the grid hU^. The S-spline of 
order 1, 

d 

(f (x) = l[_i/2,i/2]'* * l[-l/2,l/2]'* = 11'^-'^" ' 

i=i 

gives rise to the piecewise linear interpolation operator. Analogously, B- 
splines of order r G N are defined iterating this convolution r times. These 
are C""~^ (R'^) functions supported in [— r/2, r/2]'^, taking the value 1 at the 
origin. More details may be found, for instance, in 

• Band-limited sampling/reconstruction. This corresponds to the 
profile 

d 

where the cardinal sine function is defined by 

sinvrt 

sine (t) := . 

TTt 

It is easy to check that ^ (^) = 1q This profile is relevant because of 
Shannon's sampling Theorem: a function u belongs to the space 

:= |n G L2 (r'^^ : suppn C [-TT/h,TT /h)"^^ = range (^T^) 

if and only if 

In particular, such functions are determined by their values on the grid hl,'^. 

• Wavelets. Take again '■= for every /c G Z. A function 
ip E L"^ (M"^) is a wavelet provided {ipn'' : G Z'^, G Z| is an orthonor- 
mal basis of (M'^) . For more details on wavelets and the closely related 
MultiResolution Analyses, the reader may see [TIH HT)]. 

Additional examples and references (from the viewpoint of Signal The- 
ory), may be found in the survey |19j . 
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3.2 Boundedness properties 

In order to ensure that the samphng and reconstruction operators are bounded, 
we shall make the assumption (|BP|) below: 



There exist s G M and B > such that (f e H" (M) and 
r(D,)v(e):= Yl \{C + 27rky^{^ + 27Tkf <B for a.e. 



(BP) 



Lemma 3.1 Suppose (p & S' (M'^). Then the following are equivalent: 
i) if satisfies iTBPI) . 
a) There exists B > such that 



{hD,rT!;u 



(17) 



holds uniformly for h > and U ^ L?' [hTJ^^ . 
Hi) There exists B > such that 



holds uniformly for h > and u £ H 



< VB\\{hD^ 



L2 



(18) 



Moreover, whenever i), ii) or Hi) is fulfilled, the smallest constant B for 
which any of the above assertion holds is precisely {D:cYip\\l^cx,^^Qy 
Y'TOO^. To see why i) and ii) are equivalent, first observe that, given any 
9? G H'^ {W^\ the following identity holds 



T^ = {hD^)-'T!' 



(19) 



To check this, simply notice that 

tJc? (0 = /i"^ m Yl Une-'^^< = ^ {hO h'^U (hO 



hence 
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Since (-Dx)'^ (p £ L 



2 {TU>d\ 



and 



{hD^YT^lU 



L2 



L2 



it suffices to deal with the case ,5 = 0. But it is a weh-known result (see 

for instance 0113) that for ^9 G (j^d^^ 

ii) are equivalent and that 



'"</'IIloo(r'*) 



whenever is bounded. 



Statements ii) and iii) are equivalent because of the following duality 
relation: 



{hD,yT!^U,{hD,)-'v}j 



^ ' L'^{h1'i) 



which holds for every u £ H ^ 
check: 



L2(lRd) 

and [/ G (hZ'^). This is simple to 



{hD^YT!lU,{hD., 



L2 



V [/„ / (/iZ),.)^ (x) {hD^)-' u {x)dx 



(Rd)xJ?-s(]Rrf) 



Remark 3.2 For s < 0, estimate \n\ ) implies that 



{eD.,rT^U^ 



< Vb 



(20) 



a soon as h/e < 1, as it can be easily checked taking Fourier transforms. 

A sufficient condition for (|BP|) in terms of decay on ip is next given: 
Lemma 3.3 Suppose (p € CR'^) satisfies, for some e > 0, 



\{D^yip{x)\'{l + \x\f^' dx<oo 



(21) 



Then and T^^^^y^ are continuous functions. In particular, IBP|) always 
holds for such a ip. 
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Proof. It follows the lines of [151 ! Lemma II. 7. Under condition H21() . 
icy ^ £ //°'/^+^/^ (M*^); Sobolev's imbedding Theorem then ensures that 
(^)* is a continuous function and hence so is ^. The continuity of r^j^^y^^ is 
a consequence of the fact that, whenever x S (R*^) satisfies X^neZ'' Ix + 27rn)| > 
1, the expression 

r 1 

Yl ll^x(- + 27rn)||^,(^,-) 

defines an equivalent norm in H'^ (M'^) , s > 0. This actually proves that 
V sup|(0'P(e)x(C + 27rn)|2<oo. 

In particular, the series defining r^^j^^^^ is uniformly convergent and the 
claim then follows. ■ 

Condition (|2H) automatically holds for profiles <p such that 

\{D^yip{x)\ < C(l + \x\y^~^ , for every xeR'^ and some C, e > 0; 

(22) 

in particular, the hypothesis we assumed in the introduction implies 
(IbFI) for s = 0. 

Now we can prove a general result from which Proposition II .81 immedi- 
ately follows: 

Proposition 3.4 Suppose if satisfies (BP)) and we are given scales (hk), 
(efc) such that (hk/ek) is bounded. If is an h^-bounded sequence of dis- 

crete functions then the distributions m^'' [T^'-f/'^*^] are uniformly bounded 
in S' . Moreover, the limit of any weakly convergent subsequence is a positive 
measure. 

The proof this is a direct consequence of Remark 13.21 and the general 
result established in Proposition 18.11 

3.3 Bases and projections 

Below, we recall some results from Approximation Theory that will be 
needed in the sequel. These results deal with the range in (W^^ of the 
reconstruction operator T^, which we denote V^. 

The space is a Principal Shift Invariant (PSI) space. When any 
of the conditions of Lemma l3. II are satisfied, the family : n G Z*^} 

is said to form a Bessel system for V^. 
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The next Lemma clarifies how the function T(^£)^y^p characterizes further 
basis properties of the functions ip!^. 

Lemma 3.5 Let ip € S' (W^) satisfy MP\) . Then 

is an orthonormal basis of if and only if 

T(D.»(0 = 1 fora.e.C^R^. 

a) : n £ Z*^} is a Riesz basis ^ ofV^ if and only if there exist 

constants A,B > such that 

A<Ti^D^y^{0<B /or a.e. e G K'^. 

Proof. The operator {D^}" : H" (W^) {R'^) is unitary. Hence 

is an orthonormal (resp. Riesz) basis of if and only 

if l/i^'^/^ {{D^)^ if)^ : n E Z"^! is an orthonormal (resp. Riesz) basis of the 
range of T/t, >s . Thus, the Lemma needs only to be proved for profiles 
Lp ^ L? (R'^) and this is a well-known result (see, for instance, [TTl). ■ 

We shall also need the following expression for the orthogonal projection 
onto V^: 

Lemma 3.6 Let (p £ S' (M'^) satisfy IBPI) . The orthogonal projection : 
{R'^) equals = T^S^-^ {hD^Y , where, for f £ L"^ 

f £ L"^ {R'^) is defined by: 



/ (0 := 



ifTf (O/O, 



otherwise. 



This means that there exist constants A,B>Q such that 

for all U € [hU^^ . This is equivalent to the existence of a linear isomorphism R : 
V*" ^ V'^ such that | /i-'^/^i^v^;^ : n e Z"^^ forms an orthonormal basis of H" (K''). This 
property is sometimes also referred as that {'fin) „£zd form a stable frame in H" (R'*). 
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Proof. The proof of the result for s = may be found in [2] , Theorem 
2.9. We can reduce ourselves to this case by noticing that 

since, as we have seen in H19() . the range of equals that of {hD^}^^ '^(hd >V 
and {hDxY is an orthogonal mapping. Using the L^-result we obtain: 

as claimed. ■ 



4 High frequency analysis: /i ~ £ 

4.1 Reduction to the case h = e 

In this section we analyze the effect of sampling and reconstruction on 
Wigner measures at scales (e/c), of the same order of the sampling/reconstruction 
rate {hk) (i.e., such that (hk/ek) is bounded). 

First notice that it suffices to treat the case Sk = hk', the more general 
one can be obtained by a proper rescaling. This is due to the following 
identity: 

[u] {x, = {h/ef m'' [u] {x, {h/e) , 

which clearly implies: 

Lemma 4.1 Suppose hk/ek — > c > 0. Then m^'' [uk] converges in S' if and 
only if mP-'' [uk] does. Their respective limits fi^ and /x are related through: 

fi,ix,0 = c''fiix,cO. (23) 

When hk = £k, the transforms M^^ \JJ^^^^ are F-periodic in the variable 
^; hence, so are their limiting Wigner measures 

4.2 Sampling 

We start by exploring the effect of sampling on the structure of Wigner mea- 
sures. The computation of the Wigner measure at scale {hk) of a sequence of 
samples [Sj^'^Uk) is done in the following theorem; it is applicable whenever 
the hypothesis (|D]) below is fulfilled: 

essup ^ |(C + 27rn)^i^(^ + 27rn)|^ ^ as i? 

^^'3 |n|>R 
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Notice that profiles with the property H21|) immediately verify (|D]1 . 

Before stating our result, it is important to notice that the Fourier trans- 
form of a profile 99 satisfying condition ()BP|) is an element of Lf^^ {W^) . In 
particular, it is only defined modulo a set of zero Lebesgue measure. Thus, 
when dealing with pointwise properties of ^, we shall systematically assume 
that a precise representative of the class of (p has, once for all, been chosen. 

For instance, the Wigner measures of the sequences (uk) in Theorem 
14.21 below, will be assumed to satisfy conditions HMS|) and (|ND|) : 

1^ (.^ + 27rn)|^ fj-{x,(, + 27rn) , n S Z'^, are mutually singular measures. 

(MS) 

IJ. (r'^ xD^^=0, (ND) 

where, recall, stands for the set of discontinuity points of (p. These 
conditions must be understood to hold for the same representative of (p. 

Theorem 4.2 Let (hjS) he a scale and take p satisfying IBPI) and / IP)) . 

Let {uk) he a sequence in H~'^ (M*^) such that {^{hkDx)~^ u^) is bounded and 
suppose that m^'' [uk] converges to a Wigner measure fi that fulfills /TND)) . 

Then M^" [S'^'-' Uk] converges to the Wigner measure jjL^ given by: 

Z^'" 6 = E 1^ (C + 2^n)| V i + 27rn) . (24) 

Remark 4.3 i) As pointed out above, formula i28\) holds for the same pre- 
cise representative of the Fourier transform p> which was chosen in / fND|) 
and jMSJ). 

a) The necessity of hypotheses IND|) and (MS)) will be discussed in 
paragraph \4-4\ 

Hi) Condition JP)) may be replaced by the assumption that [{h^Dx)^^ u^) 
is hk-oscillatory. This will be made clear in the proof of the Theorem. 

iv) The houndedness of {hkDx)~^ Uk implies that (uk) is hk- oscillatory. 

The proof of this Theorem is postponed to the end of this section. 

The expression ()24() may be related to the concept of Wigner series 
introduced in jEHHI- Recall that given u ^ S' (M'^), the Wigner series of 
u at scale e is defined by: 
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Wg [u] (x, ^) := ^ u{x — ETm) u{x + evrn) e*"''». 

It is easy to check that Wg [u] (x, ^) = X^neZ"* '"^^ + 27rn).''' 

When is bounded in (M'^), e/j-oscillatory and possesses a Wigner 
measure at scale (e^) then the following relation holds: 

lim/ a (.,a„,-M(.,0 E'>(-.« + 2.„)*(..{), 

(25) 

for a £ S, see jS]. 

Theorem 14.21 has a simple interpretation in terms of Wigner series: the 
measure /x*^ may be obtained as the limit of the Wigner series 

Wg'" {hkDx) Uk] . 

This is due to the fact that, under any of the hypotheses (|D|) . {ip (hkDx) Uk) 
is /ifc-oscillatory. Besides, as a consequence of Proposition 18.31 the Wigner 
measure at scale (/ifc) of {'(f (h^Dx) Uk) is given by |^ (^)|^ (x, ^). The as- 
sertion then follows from l\25\i . 

As was already mentioned in the Introduction, condition (|MS|) is a re- 
striction on the support of the measure |^ (.^)|^ (x, ^). Two extremal cases 
in which it is trivially satisfied are the following: 

i) = 0, in this case ()MS|) holds independently of what fi is. 

ii) The sequence (uk) is asymptotically band-limited i.e. its Wigner 
measures at scale (hk) is concentrated on the cube Q. For those sequences, 
condition (|MS|) only involves the behavior of /i on the boundary dQ: it 
essentially expresses that the restrictions of to parallel sides of dQ do not 
overlap (i.e. are mutually singular). A sufficient condition for this is, for 
instance, 



lim sup 



^ ^0 as i? ^ oo, (26) 

{27Thkf 



where Qr := [-7r,7r - l/Rf. 

Remark 4.4 In any of the above cases, we have: 

iQ(g)/^'" {x, 6 = 1^ (01 V {x, . 

^See for the definition of tiie Wigner transform [ii]. 
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Hence, the restriction of fi^ to Mf^ x Q coincides with iJ, if and only if 
\ip (^)l^ = 1 for fi-almost every ^ € Q. 

The specific choice (p = 60 corresponds to the analysis of discretization, 
for then Sg'^u (n) = u [hn). Theorem 14.21 takes the following simple form: 



Corollary 4.5 Let (hj.) be a scale and let (uk) be a sequence in H"^ 
for some s > d/2, such that {{hkDxY Uk) is bounded. If ^ is its Wigner 
measure at scale (hk) and the measures ^ (x, ^ + 27rn) are mutually singular 
then Wigner measure fi^° corresponding to the sequence of discretizations is 
the periodization: 

In other words, /u"^" is the limit of the Wigner series w^*" [uk\ ■ 

This Corollary is particularly useful in the explicit computation of Wigner 
measures for discrete functions. As an example, consider the concentrating 
and oscillating sequences we defined in the Introduction, fk {x) = k'^l'^ p (k {x - 
and gk (x) := p{x) e'^^'^^ with p £ L'^ (R^). Using identities (0) and (jH)) we 
obtain, for (/fc) and (gi^) respectively: 

if, for instance, suppp C Q, and 

\p{x)\'dx0 ^'^+2.niO, (27) 

with no assumption on p. 
4.3 Reconstruction 

Now we deal with the reconstruction operator T^; it modifies the high- 
frequency behavior of a sequence of discrete functions in the following way: 

Theorem 4.6 Let (hk) be a scale and ([7^*=) be an h^-bounded sequence; 
take ip satisfying JEpI) . // M^'' [U^''] converges to the Wigner measure p 
which verifies JND]) then m^>' [T^^'^U^^] converges to a Wigner measure p^ 
given by: 

P^{x,0 = \^{Ofl^{x,0- (28) 
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The proof of Theorem 14.61 is based on expUcit formulas for the Fourier 
transforms of Tj^U . As we have already seen, 

T!^{0 = ^mh''U{K), (29) 

for any U L"^ (/iZ*^). The following Remark ensures that Proposition 18.31 
can be applied in the proof below. 

Remark 4.7 If £ H' {R'^) satisfies JEFl) then G {R'^; (0"). 
Proof of Theorem 14. 6L Just notice that (|29)) can be rewritten as: 

The hypotheses made on ip and /i allow us to apply Proposition 18.31 (see 
Remark 14. 7|) and conclude ■ 

Identity (|28j) expresses how the measure fj, is modulated by the profile 
the necessity of the hypothesis (|ND|) for this result is discussed in paragraph 
lO as well. 

Since fi is P-periodic in ^, formula 1)28(1 suggests that fi may be compared 
to the periodization of fi^ with respect to the variable ^. 

Corollary 4.8 Let ip, (^U^'') , and fi^ be as in Theorem \4-6\ Then the 
periodization 

:= I] + 2^n)2'^ (x, ^ + 27rn) (30) 

is a well-defineS measure, T-periodic in ^, that satisfies: 

l^^,s{xX) = T(^D^r^{i)l^{x,Cj- (31) 

In particular: 

UT{Dj;)^ip iO = 1 except for in a set of zero ^-measure then fj.^^^ = fj,. 
a) T(^[)^y^p = 1 if and only if the identity fi^ ,, = fi holds for every 
sequence {U^'^) . 

^The limit defining the sum II3UII is understood to exist for the weak convergence of 
measures in X+ (R'' x R'') . 
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Proof. Since | {^y^ ^ (^)|^ is a nonnegative continuous function, the series 
defining T(^D_^y^{^) converges absolutely for every ^ on the support of fi 
(which consists of continuity points for ^(^)). Thus, by the dominated 
convergence Theorem: 

/ a{x,(,)r^,s{OdtJ-ix,0 = Y] / a{x,(,)\{C + 2Tmy ^ + 2im)\^ dfi{x 

for every a G Cc (M'^ x R'^) . Taking now into account H28|) and the fact that 
fj, is F-periodic in ^, we find that 

/ a{x,(,)T-^,s{OdtJ.{x,^) = Y] / a{x,(,){^ + 2imf^ dfj.^{x,(, + 2im) 

and the first part of the result follows. 

Statement i) as well as the "only if" part of ii) are trivial. To obtain the 
necessity in ii), just consider sequences of discrete functions whose Wigner 
measures are of the form ^(x,^) = u (x) Snez<^ ^C°+27rn (|27|) . for 
instance). Clearly, for = fi to hold for such a measure, we must have 

r(D.)V if) = 1- ■ 

Remark 4.9 i) Because of Lemma \3.5\. if relation /x^ = ^ holds for every 
hk-hounded sequence of discrete functions then the profile (p has the property: 
^l^-d/2^h . ^ g orthonormal family in (R'^) for every h > 0. 

ii) However, the converse is not true, if ip gives rise to an orthonormal 
family then r^^^^y^ (^) = 1 holds outside a set of null Lebesgue measure. If 
fi is supported on that set, identity /i^ = /i may not hold. 

As in the preceding section, our result has an interpretation in terms of 
Wigner series. Under the conditions of Theorem 14. 6[ the measure fi^^g may 
be obtained as the limit as /c — > cxd of the functions 



hi. 



provided [{hkOxY T^'^U'^'') is /ifc-oscillatory. Note however, that this may 
not be the case for certain profiles ip (see paragraph 15. 3|) . 

In particular. Corollary 14.81 shows that the limits of vUg'' [T^^* and 
M^k ^Jjhkj coincide if we chose, for instance, (f : — -'-j_i/2 1/2)'^' 
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4.4 The necessity of the hypotheses of Theorems 14.21 and 14.61 

Formulas H28|) and (|24|) may not hold when tp is not continuous and the 
Wigner measure ^ does not vanish on the closure of the set of discontinuity 
points D^. We illustrate this with two one-dimensional examples where 

, , sin TTX 



We will chose 1q as the representative of ^ for which the counterexamples 
will be built. 

1. Necessity of condition ()ND|1 in Theorem 14.61 Take to be 
the sequence discrete function of (/iZ'^) given by their Fourier transforms: 

Then, denoting by the Wigner measure at scale h of (C/'*), 
This measure differs from /x^, which is given by: 

Remark 4.10 i) The particular choice of the representative of(p does not 
play a role. Theorem \4 b] still fails if we take as representative of the 
characteristic functions of (— 7r,7r)'^ or [— 7r,7r]'^. 

ii) In particular, this example shows that even the two projections on x 
and ^ of the measures /i and /i^ may differ. 

Hi) This also shows that the periodization in ^ o//Lt^ does not necessarily 
coincide with fi, even when = \ as is the case here. Thus the conclusion 
of Corollary \4.^ may fail when (p is not continuous. 

Our counterexample to Theorem 14.21 is essentially the same as the pre- 
vious one: 

2. Necessity of condition ()ND|> in Theorem 14.21 Define 

^{0 := l(-i,i) (^ + vr/^). 
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Then, denoting by ^ the Wigner measure at scale h of (w'^) , 



and this is different from fi'^, which is precisely: 

,2 



TT^X^ 



Finally, we investigate hypothesis (|MS|) . Now we set ip := Sq. 
3. Necessity of condition ()MS|> in Theorem 14.21 Define 

^ (0 := iq m E i^.i) - + 1) ^) • 

Clearly, as in our first example, the periodization of the Wigner measure of 
{v'') is: 

E /i (x, ^ + 27rn) = — dx ® ^ 8{2n+i)-K (0 • 

However, the sequence of discretizations {S^^v^) has the following one: 

The proof of these counterexamples easily follows from , identity (|55)) 
and Lemma 18.131 

4.5 A Poisson summation formula and proof of Theorem 14. 21 

The computation of the Fourier transform of S^^u is given by the following 
identity: 



Lemma 4.11 Let ip satisfy {BPj) and u £ H ^ (W^^. Then the Fourier 
transform of Sj^u is: 

h'' E S> in) e-'^^< = ^{K + 2^n)u U + , (32) 

ergei 
place in L\^^ 



the convergence of the first series being in Lf (W^) while the second takes 
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Proof. Begin by noticing that G (W^^ and thus 
n'^ (^) := ^ {hS, + 2TTn)u + 2tt /hn) 



is a well-defined (27r//i) Z'^-periodic L\^^ (M°') function, the series defining 
it being absolutely converg ent in L\^^ (W^) . We can compute its Fourier 
coefficients: 




Lemma l3 . 1 1 proves that S\^u is square-summable and, consequently, 

the sum being understood in the L^-sense. This is precisely formula (jc{2|) . ■ 

Remark 4.12 Identity h3l3\) may he viewed as a generalization of Poisson 
summation formula. Taking as ip the Dirac delta 6q, we obtain: 

h'^Y.^(hn)e-^'^-<=Y,nU+^-^nY 

for every u G H'^ (M*^) with s > d/2. 

Proof of Theorem 14.21 The proof will be done in two steps: 

Step 1: We first establish the result for sequences such that ^ (^) Sfc (C/^fe) 
has support in a ball B (0; R) for every k G N. We claim that the following 
formula holds: 

\n\<R+TrVd 
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This is obtained by applying the inverse Fourier transform at both sides of 
identity ()32|). and remarking that only summands satisfying \n\ < R + Ti\fd 
must be considered because of the condition on the support of 'ipv^ {■ /hk). 
The Wigner measures of the functions 

e-^^'^--'^^^{hkD^)uk {x) 

are precisely (cf. Proposition 18.31 and Remark 14.71) : 

|^(^ + 27rn)|V(2;,^ + 27rn) . 

By hypothesis, they are mutually singular so, by Lemma f8.131 we deduce 

^So '^'P 



is given by 



that the measure /j^ obtained as the limit of m^'' 
(El- 

Step 2: We prove the result in the general case by taking advantage of 
hypothesis (|D|) . Let x ^ CU.'^) be a cut-off function identically equal to 



one in the unit ball B (0; 1) . Denote by S^''j^Uk the truncation given by: 



R 

^ + 27rn 
hk 



Then, by the first step we have just proved, M''' 



converges to 



^ + 27rn 
R 



^ (x,^ + 27rn) 



We claim that (|D]) implies the following: 

2 

L2(hkZ<i) 



lim sup 

fc— >oo 







as ^ oo. 



It is sufficient to realize that 

for ipfi := xi'/R)'?- The norm of is precisely (cf. Lemma IXT|) 



essup > 

«6Q ^ 



+ 2'nnY + 2im) x 



^ + 27rn 
R 



(33) 



(34) 
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which tends to zero as R ^ 0. 

Lemma [8. 121 then ensures that /i^ weakly converge to /i"^. Identity 
means that 

/ a{x,i)diil{x,i)= ! y^a{x,i + 2^Tn)\^{i)\'\x{i/R)\'d^,{x,i) 

for every test function a € S. Passing to limits as ii ^ cxd in the above 
identity we obtain the claimed result. 

Notice that the same argument may be applied if, instead of condition 
(|D|). we have that {{hkDx)~^ Uk) is /i/j-oscillatory. This is because may 
be estimated from above by 

2 



lim sup 



{hkDx) 



1-X 



R 



Uk 







as ^ oo 



L2 



because of Lemma l3. II and the /ifc-oscillation hypothesis. 



5 Computation of defect measures 



5.1 Relations between defect and Wigner measures in the 
discrete setting 

In this paragraph, we establish the analog of Proposition II .21 in the discrete 
setting. In particular, we present conditions that ensure that the projection 
on the x-component of a Wigner measure may be obtained as the limit of 
quadratic densities of the type: 

2 



(x) := h'^Y.K 



'hn 



{x). 



Proposition 5.1 Let (hk) be a scale and be an hk-bounded sequence. 

Suppose that [M^^ [U^^]) converges to as k ^ oo. Then, for every (f) G 

Cc 



dfi{x,C) = lim (hk)'^ 

fe— >oo ^ — ' 



kn 



Tjhk 



(35) 



If (sk) is a scale such that hk <C £k o,nd the transforms M^* converge 
to fi, then / l^5|) holds provided (U^'=) is s^- oscillatory, i.e.: 



lim sup (hk)" 

fc— >oo 



Q\B(0;/ife/efei?) 



Uhk (^) 



dS^ ^ as R ^ oo. 



(36) 
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In view of Proposition 15.11 one could think that Wigner measures at 
scales coarser than h/^ are unnecessary. However, as the next result shows, 
if is efc-oscillatory for such a scale then the Wigner measure at scale 

(/ifc) does not give any information about the oscillation effects. 



Proposition 5.2 Let (hk) and (e^) be scales such that £k- For every 

Ek-oscillatory, hk-bounded sequence such that M^'' [f7'^'=] ^ fj, as k ^ 

oo we have 

11 (X, ^) = (x) ® ^ d2nk (C) 



where v is the weak limit in M. 



') ofth 



e measures 



The Wigner measure also gathers the information on the densities | J-'^'= U^'^ (0 1 ^ i 
indeed, these converge to the projection on ^ of the Wigner measure provided 
that no energy is lost at infinity. 

Proposition 5.3 Let (hk) and (e^) be scales such that hk/sk is bounded. 
Suppose that {U^*') is compact at infinity: 



limsup(/ifc)'^ ^ |c/^ 

\hkn\>R 



k^oo 



as R 



oo, 



an 



d that W"^ [C/^fc] ^ /i as k ^ CO. Then 



I 



lim 



di 



(37) 



for every £ Cc 



The proof of Propositions 15.11 and 15.31 requires the following preliminary 
result, which explains how the transform of a discrete function 

can be localized: 



Lemma 5.4 Let G (/iZ'^) and (f, (/) £ C~ (M"') . Then for every 
aeS (W^ X R'^) the following holds: 



lim 



Jjh,k 



S'xS 



hk J {2TT f 



0. 
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Proof. First remark that, as a consequence of relation (|14p and Lemma 
18.51 we have 



Mm 

fc^oo 



S'xS 



S'xS 







(38) 

for every test function -0 G {W^') such that -0 (x) = 1 for x G supp (j). 
Now, (|521 i) and (|14jl together with Plancherel's formula for the discrete 
Fourier transform yield: 



S'xS 



hkTThk 



S'xS 



\2d 



^ 2 

(27r) 



and the result follows. 



Proof of Proposition 15.11 Identity (jSHJ) in the case = is a 
direct consequence of the identity 



/ 

JQ 



M''[U]{x,Od^ = E'' (x) 
and that, due to the F-periodicity in ^ of M^'' and fi, one has 



lim 

fc— >oo 



'xQ 



(j){x)dn (x,0 



for every cp e C° 

Next we analyze the case hk/ek 0. Given functions 0, x e (K'^), 
using Lemma 15.41 and periodization in the variable ^ we find: 



Choose a function x £ C'^ iW^) such that 



X 



(39) 



X (0 = 1 for lei < 1, 

X (0 = for lel > 2, 

< X (0 < 1 for C G IK'^, 
and set Xr (0 •= X iC/R) for every R > 0. With such a test function and 
/ifc/efc < vr/i? we have Xi? (C + 2vrn)^ = Xr [j^C) < 1 for every ^ S Q. 
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Then, taking this into account in H39|) and using Plancherel's formula, the 
following is obtained: 



lim 

fc— >oo 



(27r)^ 



where 



M{R) := limsup(/ifc)'^ 

fc— +00 



1 - Xr ( 1^ + 27rn) ) ) 0C/^^ (0 



< M (R) 



Identity is obtained by letting R tend to co, noticing that implies 
that M {R) ^ Q as R ^ oo. m 



Proof of Proposition [5T2l Since 



km 



is e/c-oscillatory, we have 



lim sup 

fe-^oo JQ\B{Q;S) 



for every 6 > and cf) G (M'^) . Using Lemma 15.41 below, we obtain, for 
every ^gC-(Q\{0}), 



k—>oo 



0= lim / ifiC) ^U^k{^) 



R^xQ 



In particular, is concentrated on the set x {0}. Since /U (• x Q) = i/ (x) 
by Proposition 15.11 we find that, because of the periodicity, fi (M'^ x •) = 
u (M'^) YykeZ'^ "^aTrfc (C); and this restricts to be equal to u (g) Y^kez^ ^2-Kk- ■ 



Proof of Proposition ISTSl Since the densities are uni- 

formly bounded in (W^^ (and consequently in M. {W^') ) it suffices to prove 
the result for test functions V £ C'^ (M'*). Let x and Xr be defined as in 
the proof of Proposition 15.11 Because of Lemma 15.41 the following holds for 
every G 5 (W') : 



\Xr{x)\ ij{Odfi{x,0 = lim 

fe— »oo 



d^. 



Since \xr ix)\ 1 as R ^ oo for every x G M'^ we only have to show that 



lim 

R—->oo 



\XR{x)\^'^lJiOd^l{x,0= lim / V (6 J'^'U^" (0 

k--*oo 



d^. 
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This appears as a consequence of the identity 



which imphes 



Hm sup 

fc^oo 



^^^^ r-.u^'^ii) - j'^'^xrU'^HO 



di 



< Cw, hm sup 



u 



Since the U^'' are compact at infinity, the second term in the above estimate 
tends to zero as R tends to infinity and thus: 



hm sup / V' (0 

fc— >oo 



as 



One easily deduces from this that the measures \j^^''U^'' converge 
in (R'^) to the measure J^a A* {dx, •) as claimed. ■ 



5.2 Defect measures of reconstructed sequences 



Let (f/'^*) be /i^-bounded and ip G (M*^) some profile satisfying (|BP() . 
As a consequence of Lemma 13.11 the sequence of densities 



is uniformly bounded in (M'^) . Hence, Helly's compactness Theorem en- 
sures that, extracting a subsequence if necessary, there exists a measure 
£ M+ {R'^) such that 



lim 



ix) 



dx 



R<*xQ 



4> (x) dv^ (x) 



The main issue addressed in this section is that of clarifying how u^p depends 
of the sequence and the profile (p. We shall see that a formula relating 

and the limit of E^'' [C/'^'=] does not exist in general. However such a 
formula may be established in terms of the Wigner measure of (f/'**). 
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Suppose that M'^'' [U'^'=] converges to fi. Then, Theorem 14.61 may be 
appHed to obtain that, provided (M^ x D^) = 0, one has 



m 



In general, we are only able to ensure (see Proposition 1.7 in jH]): 



1^^ [x) > 



but equality holds whenever (^{hf^Dx)^ T^''U^'=^ is /i^-oscillatory. However, 
this is not immediate since there exist profiles ip for which {^{hkDxY T^'^U^'') 
may fail to be /ifc-oscillatory for some (an example is provided at the 

end of this section). Nevertheless, 

Proposition 5.5 ((/i^Z)^)* T^'^f/'^*) is hf^-oscillatory whenever JP)) holds. 

This immediately follows from: 
Lemma 5.6 (^{h^Dx}^ T^'^U^''^ is hk-oscillatory if and only if 



k—>oo 



d^ ^ as R ^ oo, 



where 



\n\>R 



Proof. Start by noticing that: 



Per iodizing in we get: 



d^ 



Q 



(0 



[//ifc (^) 



d^ < 



d^. 



d^< I aj-^'^(e) c^(e) 



and the claim follows. 



Hence, /ifc-oscillation is obtained if decays at infinity at an uniform 
rate. For more general ip, it is still possible to obtain sufficient conditions; 
however, these depend on the particular sequence of discrete functions to be 
reconstructed. 
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Proposition 5.7 Suppose 

i) fi (r'^ X D 



0, 



ii) (U^'') is compact at infinity. 



(40) 



Then {{hkD,j,Y T^'^U^'') is hk-oscillatory. 



Proof. For the sake of simplicity, we prove the result for s = 0; the proof 
in the general case being identical. Taking into account the periodicity of 
the densities involved, Proposition 15.31 ensures that 



lim (hk) 

A;— >oo 



V' iO u^^ (0 



and the claim follows. Since /x CR.'^ x D-r 



(41) 



0, 



for every Cc 

necessarily /x ( R"' x D^r ) is null for every R > 0. Prom classical results on 



weak convergence of measures, one deduces that relation ()41() also holds for 
ijj = a^. Hence, by the dominated convergence Theorem, 



lim lim 

R-^oo fc— >oo 



(hk)' [ (0 
Jo 



Uhk (^) 



lim 



a« (0(i^(x,0 = 0, 



and the result follows. 



What we have proved so, combined with Prop osition 1 1 . 2 1 far is gathered 
in the next proposition. 

Proposition 5.8 Suppose at least one of IP]) or {^ifjj is satisfied and that 
fi {R'^ X D^) = 0. Then 



lim 

fc^oo 



dx 



for every (j) £ Cc 



<Pix)m'^iordf,ix,o 

(42) 



As anticipated above, H42|) shows that the knowledge of the weak limit 
of the measures E^'' ^U^''~\ and the profile (p are not enough, in general, to 

reconstruct the weak limit of the densities | (/i^-Da;)* T^* |^ dx. However, 
when the sequence of discrete functions under consideration is efc-oscillatory 
for some scale coarser than the reconstruction step /i^, there does exist a 
formula that relates both limits: 
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Corollary 5.9 Let {U^'=^he an hk-hounded, e^- oscillatory sequence such 
that [E^k [t/'^ft]) weakly converges to a measure v. Suppose moreover that 
ip is continuous at T and that any of ID)) or i4(J{ ) is satisfied. Then the 
densities \{hkD^yT^<-U^>^\ weakly converge to the measure 



Proof. Using Proposition we find that any Wigner measure at scale 
hk of equals 

H (x, =v{x)® ^ 52-nn (C) • 

Since (M'^ x D^^ = 0, Proposition 15.81 is applicable and gives: 

{hkD^Y r^'=[/'^'= I' dx - T^a^y^ (0) V (x) as ^ oo, 
as claimed. ■ 

Remark that condition (|4Ul i) reduces in this setting to the requirement 
that T(^j^_^y^ is continuous at ^ = 0. 

5.3 A counterexample to /i-oscillation 

Here we exhibit a function tp ^ LP' (M) satisfying ()BP|) with p is continuous 
but 

Ik^lL-CQ) = ^ ^^^"^y i? > 0. 

With such a profile, we show that there exist a sequence of discrete functions 
{U^) such that (T^U^^ is not /i-oscillatory 

To construct ip, define t„ := e"" for n = 0,1,2... and let ip^ be the 
piecewise linear function given for n > 1 by 



t — tn+1 
t — tn-1 



if t G (t„+i,t„) , 

if t G (t„,t„_i) , 



in tn~l 

otherwise. 
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Clearly X^^i V'n (^) = 1 for t G (0,ti), and the sum vanishes for t < 0. 
Defining 



n=l 



we obtain tp € (M'^), ^eC {R'^) and (^) = ^^^^ (6- 
Moreover 

■ 1 ifee (o,x„+i), 



1 for every > 0. 
If we chose discrete functions G (/iZ) such that 



Thus 



(0 = h-' 1(0,/,) (e + 27rn) 



ne2 



then for the if above constructed we obtain 



hm / (^) h {(,) ' = 1 for every R>0. 



Q 



This proves that (T^U^^ is not /i-oscillatory. 



6 High frequency analysis: /i <C e 

Here we shall investigate the structure of Wigner measures at scales (sk) 
asymptotically coarser than the sampling/reconstruction rate (/ifc)- 

In the next two Theorems, we suppose that (p satisfies (|BP|) and (p is 
continuous in a neighborhood o/ ^ = 0. Moreover, (hk) and (sk) will be 
scales such that hk <C Ek- 

Theorem 6.1 Suppose that (U^''^ is h^-hounded and M^'' [f7'^'=] converges 
to the Wigner measure fi. Then m'^'= [T^* converges to a measure fx^ 
given by: 

fi^{x,o = \^m^t^{^,o- (43) 

The proof of this result is completely analogous to that of Theorem 14.61 
Concerning the sampling operators, the situation is much similar: 
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Theorem 6.2 Let (uk) be a sequence in (M"^) such that {{hkD^y Uk) 
is bounded in (W^) and ek-oscillatory. 



i) Then [S^'' Uk) is Ek-oscillatory. 

ii) Suppose moreover that m^^ [uk\ converges to a Wigner measure fi. 
Then M^'^ [S^>^ Uk\ converges to the Wigner measure /x*^ given by: 



/x'^ = 1^(0)1 V- 



(44) 



Proof. To prove the first part of the Theorem, begin by noticing that, 
by the Cauchy-Schwarz inequality and Lemma [4.111 for almost every ^ G M'^: 



S^'^Uk {Cl 



< 



{hk) 



2d 



—d^^{i + 2^n)uk\-^^ 

{i + 2Tm)-'vrk 



Thus 



So 



{hk 



i + 27rn 
hk 



Sip'^uk {0 



Q\B{0;R/et:) 
< llr, 



Q\B{0;ht:/6kR) 

|(/ifee + 27rn)-'u^(e + 27rn)|'de 



2 



\{hkO^'uk{0\ dc, 



lR''\B{0;R/sk) 

and this clearly proves that {S^'^Uk) is efc-oscillating as soon as {{hkDx)~^ Uk) 



IS. 



The proof of identity (|44|) is essentially identical to that of Theorem 
14.21 A completely analogous argument to that used in the step 2 of that 
proof allows us to consider only sequences such that (f {hk/£k')uk {'/^k) is 
supported in a ball B{0;R). This hypothesis together with Lemma |4.11l 
implies that, for hk/sk small enough. 



{hk)'^ S^'^Uk ( — ^ ) = P ( — ^ ]uk{i/ek), 



hk, 



hk, 



that is, only one summand is involved. Then the result follows from Propo- 
sition exactly as in the proof of Theorem 14.21 ■ 

We conclude with a simple remark: 
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Corollary 6.3 Under the assumptions and notations of Theorems \6. 1\ anc 



i) If ip has zero mean (i.e. ip{0) = 0) then the Wigner measure at 
scale (sfc) of any sequence (T^'=[/'*'=) or (^Sj^'^Uj^^ vanishes identically. In 
particular, this is the case if is a wavelet.^ 

ii) fj,^ = fi"^ = fj, always holds for profiles such that |^ (0)| = 1. 

7 Wigner measures of Sampled/Reconstructed se- 
quences 

Now we are able to describe Wigner measures of sequences of the form 
T^5^n. In its full generality, our result requires several compatibility hy- 
pothesis, that we describe below. First of all, 

i) V and ip satisfy (|BP|) with exponents s' and s respectively. 

(45) 

ii) if satisifes (|D|) . 

The admissible sequences will be assumed to be such that: 

Uk G H-" (r'^^ and {{hkD^y Uk) is bounded in (r^^^ , (46) 

and their Wigner measures must satisfy the following compatibility condi- 
tions for some precise representatives of ip and ^: 

i) fi fulfills (INHll . 

ii) / l^(e + 2™)|i^(0|'dM(x,0 = 0, neZ'^. (47) 



iii) /i satisfies ()MS|1 . 
Combining Theorems 14.61 and 14.21 we obtain: 

Theorem 7.1 Let tp and ip he functions satisfying ( |^5| ); let (hk) be a scale 
and {uk) be a sequence satisfying Suppose moreover that m^'= [uk] 

converges to a Wigner measure that satisfies (^^. 



Then m^^ 



rpH qhk 
-'V "^f 



converges to the measure ^^^^ given by: 
[ a{x,C)dn^^^{x,0 = [ V a(x,^ + 27rn) ^(C + 27rn) \^p{(,)\^ dfi{x,0 



(48) 



'See, for instance, ITDi . Proposition 2.1. 
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for every a £ Cc 



Proof. Hypothesis (|47l iil expresses that the closure of the set of dis- 
continuity points of V) is a nuh set for the Wigner measure of S^^Uk, 
X^fcGZ'* 1^ + 27rn)|^ /i (x, ^ + 27rn). Hence, Theorem 14.61 is apphcable and 



we conclude that the distributions m^*' 



converge to the measure 



Since V iO is integrable with respect to the finite measure \ip\ /u (this is 

again due to (^Tjii)), its periodization is integrable as well and formula 
follows. ■ 



Remark 7.2 i) When ip and verify 1121]} hypotheses \45\ ), \4T\i) o,nd 
i4T\ ii) are immediately satisfied. 

a) i45\ ii) may be replaced by the requirement that {{hkDx)~^ Uk) is h^- 
oscillatory. 

From formula (|48j) one sees at once that, taking ip = ip = 5q one has 
that fi^^^ is the periodization in ^ of the Wigner measure /i. Hence, /i^ ,^ 
coincides with the limit of the Wigner series corresponding to (n^). 

When ip and |^|^^ vanish off Q it is easy to check that formula ()48|) 
takes the simple form: 

^^^,^(,{x,0= 1^(61 i^ix,0- 

It is also clear that, as soon as |^ (C)l^ (^) is not null outside Q, the 
measures fj,^^^ and fi will in general differ. 

Concerning defect measures, combining ProDosition l5.8l and the previous 
theorem, we obtain: 

Theorem 7.3 Under the notations of Theorem \7.1\ the following holds: if 

s' h h ^ 

(hkDx)^ T^'^S^p'^Uk dx weakly converges to a measure u^p^^ 
and ijj verifies IP]) then: 

i^p,^{x)= I Yl |(C + 2^n)^'^(e + 27rn)|'|i?(e)|V(^,de)- (49) 
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Remark 7.4 i) The conclusion of the Theorem still holds if condition "ijj 
satisfies IP])" is replaced by \4(^ - 

a) Theorems \l.'A and \1.5\ follow immediately from Theorems \7.1\ and \7.'J\ 

With formula (|49|) at our disposal, we are now able to answer, in a quite 
general way, the questions A-D addressed in the introduction. Of course, the 
answer to A is negative, since, in general, /j, is not trivial in its ^ component; 
concerning the problem of filtering, we immediately get the necessary and 
sufficient condition: 

I 2 



Analogously, c^,^('?) = 1 for /i-a.e. S R"^ characterizes the profiles that 
give v^^ = V. To answer D, we must, of course, assume that p and r , ,5/7 

are continuous (which, as we know, is the case if ()21|) holds). In that case, 
we have and equality v^p^^ = v for every admissible sequence if and only if 

The sampling profile cannot be an (R*^) function, since |p|^ is neces- 
sarily periodic. When Lp = and il) generates an orthonormal basis in the 
sense of Lemma 13.51 we always have Vip^^i) = v. If merely generates a Riesz 
basis, Av < v^p^^ < Bv holds instead. 

The above results may be used to compute Wigner measures of the or- 
thogonal projections P^''uk of a given sequence {uk) on the shift-invariant 

space defined by the range of r;^ As we have seen in Lemma 13.61 may 
be written as the composition of with S"^ {HDxY for a sampling profile 

(p := {DxY ip- Hence, Theorem 17. II gives: 

Corollary 7.5 For ip satisfying \21\) and (n^) such that (|^6p and / TlVrSj) 

holds, the defect measures of the sequence (Plp'^Ukj given by: 



^^^^^ 



where 1^ (^) denotes the characteristic function of the set of ^ such 
that Ti^D^Y^ (^) / 0. 
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In particular, when ip gives rise to an orthonormal family, we obtain the 
simple formula (cf. Lemma l3.5() : 



To conclude, we shall see how the above results may be refined when 
the sequence {{kkD^)''^ Uk) is assumed to be efc-oscillatory at some scale 
hk £k- The assumptions of ip and are weaker: 

i) il) and (p satisfy (jBPp with exponents s' and s respectively. 

ii) ip, ^ are continuous in a neighborhood of ^ = 0. (50) 

iii) T I .s' is continuous at ^ = 0. 

Theorems 16. 11 and Corollary 15.91 give then: 

Theorem 7.6 Let ip and (p he functions satisfying \5IJ\) : let {hk), (efc) he 
scales with hk ^ £k o.'f^d let {uk) he a sequence such that \4b\j holds and 
(^{hkDx)~'^ Uk) is Ek-oscillatory. Suppose moreover that m^'= [uk] converges 
to a Wigner measure fi. 



Then m^'' 



'Uk 



converges to the measure fj,^^^ given hy: 



Moreover, if 
then: 



{hkDyTl;'=s!^^ 



'Uk 



dx weakly converges to a measure 



i^</3,v(^)= |(27rn)'' i){2TTn) |^(0)|^z^(x) 



where v is the wGCik liiTiit of the dcTisitics ^{h]^Dx 

) ^ Uk\ dx. 

Hence, when a sequence possesses a characteristic oscillation scale {sk) 
(that is the meaning of the e^-oscillation condition), choosing a sampling/ 
reconstruction rate {hk) asymptotically finer than {sk) allows to completely 
capture its oscillation/concentration behavior (modulo a constant that only 
depends on -0 and ip). 

Filtering in that case can only be achieved by means of a sampling profile 
ip with zero mean {p> (0) = 0) or a reconstruction profile such that vanishes 

at r. 
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8 Tools from the theory of Wigner measures 

The main tools from the theory of Wigner measures used in this article are 
Propositions 18. ll and l8.3l below. The first of these is an extension of Theorem 
ll.ll to bounded sequences in Sobolev spaces: 

Proposition 8.1 Let {sk) be a scale and (uk) be a sequence of functions in 
j^~s (^d^^ some s > satisfying: 

\\{£kDx)~^ Uk\\^2(^d^ 0''^^ uniformly bounded in k. (51) 

Then the sequence of distributions [m^^ V^k]) is uniformly bounded in 
S' . Moreover, any of its weakly converging subsequences tends to a positive 
measure. 

As we have done so far, a measure ^ G (R'^ x E"^) will be called 
the Wigner measure at scale (e^) of a sequence (n^) (satisfying the 
hypotheses of Proposition 18. 1|) provided m"^* [uk] ^ in S' as k ^ oo. 

Remark 8.2 i) When s > 0, condition l\51]) is stronger than just requiring 
that (uk) is bounded in H^^ (M'^) . 

a) Let (/ifc) be a scale such that hj. <^ £k- If \\{hkDx)~'^ Uj.\\^2^-^d^ ^ C 
for every /c G N then \\{£kDx)^'^ Uk\\^2(^^d'j is uniformly bounded as well. 
Hi) The same result holds if [■] is replaced by the Wigner transform 

The second main result of this section a localization formula for Wigner 
measures which is used several times in this article: 

Proposition 8.3 Let (ek), [hk) be scales and let (uk) be a sequence in 
^-s (i^cf^^ s > 0, satisfying k51\). Suppose that (j) is a Borel function such 
that (j) G L°° (R'^;(^)'^) for some r G M. If m^'' [uk] converges to n then 
nf^ [(/> {hkDx) Uk\ converges to a Wigner measure /x^ which has the follow- 
ing properties: 

i) If hk = £k cind (R'^ x D^^ = 0, D^f, being the set of points where (p 
is not continuous, then 

H (^>6 = l'A(0lV(^,0- 

a) If hk ^ £k o-n-d (p is continuous in a neighborhood o/C = then 

M^ = l</'(o)lV- 
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When applied to := (0*, this result gives: 



Remark 8.4 Let (e^), (/ifc) and (uk) be as in Proposition Suppose 
m^'' [ufc] converges to fi. Then m^*^ [(/ifc-Ca;)"* ^fc] converges to the measure 
given by: 



Hs = fJ; if hk<^ Sk- 

in particular (cf. Theorem \l.l\) . n (resp. fi) is a finite measure when 

hk = (resp. hk <. Sk). 

For the convenience of the reader, we give detailed proofs of both results; 
they follow the ideas present in the existing literature on the subject ( |6lll3| 
ISlini)- Prop osition 18 ■ 1 1 will be proved in paragraph 18.21 We shall essentially 
show that truncation of the high frequencies of a sequence satisfying (|5H) 
implies ,^-variable localization of the corresponding [•]. Then we conclude 
by applying Theorem II. II to the localized sequence. 

Proposition 18.31 is proved in paragraph 18.31 to conclude this section, we 
describe two results useful for the computation of Wigner measures (Lemmas 
lO^ and imil . 

8.1 First properties of [u] 

We begin by discussing three alternative ways of computing [u] that may 
be used when u is merely a tempered distribution. First remark that, given 
a u ^ S' (M*^) , it makes sense to consider the distribution [u] given by 
((21), since the Fourier transform of u is well-defined. Actually [u] £ S' . 

1. The action of [u] on a test function a G 5 is given by any of the 
formulas (see 0): 



(m^ [u] , a) c, 



S'xS 



{u,a{x,eD^)u)^,^^a^^s(Rd^, (i) 
(x, - — -] u {p) u {x)dpdx. (ii) 



(52) 

where a{x,eDx) is the semiclassical pseudodifferential operator of 

symbol a: 

a (x, sD^) u{x)= f a (x, eS) u (C) e'"'^-^, (53) 
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and the kernel ka ix,p) is the inverse Fourier transform of a with respect to 



kaix,p) := / a{x,0^ 



(2vr) 



Formula (|52l i) makes sense because the operator a{x,eDx) maps contin- 
uously S' (M"^) into S (M'^) whenever a G 5 (see, for instance, CHI)- The 
integral in (|521 ii') must, of course, be understood in distributional sense. 

2. The distribution [u] may be computed through the rescaled Fourier 
transform 



(2vre) 



d/2 



using the identity: 



m'^ [u] (x, ^) = [J^^u] (^, —x) 



(54) 



(55) 



This follows from a direct computation from the definition (j2}. 
3. Now we present two localization formulas: 

Lemma 8.5 Let u e S' (M°') , E C°° {W^; {xY) for some r eRandaeS. 
Then there exists r",r2 G 5 such that: 

(m^ ,0)5,^5 = (l'A(a;)|^m= N ,a) +e(m=['«],rf)5,x5' 
(m^ [0(/iD^)n] ,0)5,^5 = - 



e 



m \u\ , a 



S'xS 



e \ ' ' ^ /S'xS 



Moreover, the test functions rf, are uniformly bounded in S for < a < 
1. 



This holds as a consequence of standard results on symbolic calculus 
for semiclassical pseudodifferential operators; see for instance Remark 
that Proposition 18. HI is not a consequence of this result, since the multiplier 
(j) (hDx) there may have a non-smooth symbol. 

8.2 Boundedness of the transforms [u] 

The next lemmas are used to establish the boundedness in S' of the sequence 
{m^'' [uk]) provided (uk) satisfies the hypotheses of Proposition 18.11 
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Lemma 8.6 For every u ^ L"^ (M'^; (x)*") and a ^ S the following estimate 
holds: 



|(m^ [n] ,0)5,^^1 < ||n|| 2/ / sup ka{x,p) {x - ep) {x) 

Proof. Use formula ()52lii) to write 

{rrf [u] ,a) gty^g = I / ka{x,p)u{x — ep)u{x)dpdx 



dp, 



noticing that this integral makes sense as ka £ S. Multiply and divide the 
integrand above by {x — epf^'^ {^Y^"^ to obtain, by Holder's inequality, 



dxdp, 



where we have set Ur (x) := {xY^'^ u (x). The conclusion follows from another 
application of Holder's inequality. ■ 

If tt E H-" {W^) then J^^n G (w^-^ • Clearly, 



\{eDx) '''"||^2(]8d) = ll^^^^n"^ 



(56) 



Thus, taking identity into account, we obtain using the preceding 
lemma: 



{rrf [u],a)g,^A < \\{eD^) 



L2 



/ sup|S(g,e)(e + eg)'(0^ 



dq 



{2^r 

(57) 

where a {q, ^) denotes the Fourier transform in x of the function a (x, ^). 
Lemma 8.7 For every s > there exists a constant Cg d > such that 



1 1 2 



holds for every u £ H {W^\ and every a £ S 



sup 



a{q,i){iY' {eqYdq, 
(58) 



Proof. This is obtained through the simple inequality {£, + q)^ < C^^d {(,)^ {q)' 
which holds when s > 0. ■ 

Notice that whenever a £ S, the integrals J^a sup^g^d a {q, ^) (0^* (^Q')* dq 
are uniformly bounded for < e < 1. Consequently, 
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Corollary 8.8 Let (e^) and (uk) satisfy the hypotheses of Provosition \K7[ 
Then the sequence [rrf^ [^k\) "is bounded in S' . 

Estimate ()58p gives immediately the following: 

Remark 8.9 Lemma \8. 7| shows that [u] acts continuously on test func- 
tions a in the closure of S for the norm: 

[al := / sup a{q,C) (O'' {q)' dq < oo. (59) 



This closure contains the space 

:= {(D^)^ (O^s ^^^^ l^^d ^ jgd^ . [^j^ ^ ^1 ^ (gQ) 

Remark 8.10 Consequently, if {uk) is as in Provosition \8.1\ and {rrf' [uk\) 
converges weakly in S' then {m^^ [u\ ,a) converges as well for every a £T,^. 

Proof of Proposition l8.1L The boundedness of the sequence (m^'' [ufc]) 

was proved in Corollary 18.81 Suppose now that the distributions m^'-' [uk] 
weakly converge to some n G S' . We next show by means of a localiza- 
tion argument that is a positive distribution and thus, due to Schwartz's 
Theorem, a positive Radon measure. 
Take (p e S (ufj ; Lemma lO gives 

lim (m^'= [(j) {£kDx) Uk] , a)^,^^ = a{x, ^ \(p (01^ dfi {x, 

for every a £ S. Since {4> {ekDx) Uk) is a bounded sequence in 
Theorem 1 1 . 1 1 ensures that \4>{S,)\'^ fi is a positive Radon measure (and hence 
a positive distribution). But 4> G S (^fj is arbitrary, so itself is positive 
and we obtain the desired result. ■ 

Notice that a very similar proof would give a version of Proposition 18.11 
in the context of weighted spaces (M'^; (x)''). 

8.3 Proof of Proposition lOl 

The key ingredient in the proof of the Proposition is the following auxiliary 
result: 
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Lemma 8.11 Under the assumptions of Provosition W^ and for every a G 
S, if any of the following conditions hold: 

i) hk = £k ^iT'd a vanishes on the set of discontinuity points of (j). 

ii) hk <^ £k o-nd (p is continuous at = 0. 
Then 



lim 

fc— >oo 



nf - [uk] , a 



0. 



(61) 



Proof. Take a € 5 and set (0 •= {^k/ ^ki)- Prom relations ((^, 
H52li) and we obtain: 



m^'= [<i){hkD^)uk] 
where 



hk^ 
£k 



nf^ [uk] , a 



< Mk (a) ll(efe-Dx) Uk\ 



L2(Kd) 



Mk (a) := / sup |a (g, ^k (C) [^k + euq) - ^k (C)] + euqY {iY 



dq 

(2^' 
(62) 

recall that o (g, ^) stands for the Fourier transform of a {x, ^) in x. 

We now must prove that Mk (a) — > as A; — > oo. This will be done by 
first checking that for test functions a belonging to the smaller class: 



V 



a £ S -.a £ C, 



OO / Tad ^ TTDd 



)}■ 



Take R > such that suppa is contained in B (0; R) x B (0; R). 
When A: G N is sufficiently large, ^ 1 and 

— {C + £kq)^B{0;2Rsuphk/ek) for every g, ^ G 5 (0; i?) . (63) 

Suppose now that i) holds. If denotes the set of points where (f) is 
continuous, then <I>fc = is uniformly continuous over C^j^ Cl B (0; R) and, 
consequently, 

sup lc^(6l'/'(C + efc9)-</'(6l ^0 as/c^oo 

g,5eB{0;i?) 



because of (jHSJ- 



50 



On the other hand, when hk/sk ^ and cj) is continuous at ^ = 0, again 
as a consequence of (|63() . 



sup 

C,<?6i?(0;i?,) 



<2 sup |0(O-</'(O)| ^0 

5eB(0;2/ifc/£feiJ) 



Thus, in either case, 

sup \a (q, e) ^k (6 [^fc + e^g) - ^fe (C)] + efc9>' (01 ^0 as A: ^ oo 



for every q G M.'^. Lebesgue's dominated convergence Theorem gives the 
convergence to zero of the integrals ()62() . The density of I? in 5 concludes 
the proof of the Lemma. ■ 

Proof of Proposition 18. 3L To prove i) and ii) it only needs to be 
checked that, for any a £ {R'^ x W^) (if Sk = h^, we further require that 
a\^dy^jj^ = 0), the functions |(/> (/ifc/efc^)^ a (x, ^) belong to the class T,^. If 
so, then 



lim 



S'xS 



holds with c := lim/ifc/e^., because of Remark I8.1UI The conclusion would 
then follow from identity ()61() . 

First, notice that |(/> (/ifc/efc-)!^ a are compactly supported and infinitely 
differentiable in x. When Sk = hk we must verify that |(/>|^a G which is 
clearly the case if al^dy^jj- = 0, for then |(/)pa is continuous in ^. 

On the other hand, if hk <^ and <j) is merely continuous in a ball 
B{0;5) then, for k large enough, suppa C B{0;hk/skS) and consequently, 
(p{hk/£k') is continuous on suppa. ■ 



8.4 Additional properties. 

The next approximation result is sometimes useful in the computation of 
Wigner measures: 

Lemma 8.12 Let (uk) and {u^) be sequences in H-" iW^), s > 0, satisfy- 
ing \51\) with the same hound and 

\i-ms\xp\\{£kDxy {uk - Uk)\\^2(^d) ^0 as iV ^ oo. 
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Suppose thatm^'' [uk] andm^'^ [u^] converge respectively to ^ and ji^. Then 
fiN^ 1^ in M+ (^'^ X M"') as N ^oo. 
Proof. This is a simple consequence of the identity: 

1'^'= [uk] - m^'^ [u^] , a)^,^^ = /u^,a (x, SkD^) (n^ - n^) ) + 



m 



+ {[Uk- u'-^), a (x, SkD^) Uk 



S'xS 
S'xS 



This gives an estimate: 



{m'^iukj-m"^ [u^] ,a)^,^^ <C\\{ekD^) ' K - <) ; 
taking hmits as /c ^ oo we obtain: 



a{x,^) (dfi-d^^) 



< Chmsup ll(efc-Dx)'' {uk - Uk)\\L2( 

k^oo 



and the result follows, since the measures equibounded. 
We conclude this section with an almost orthogonality result: 



Lemma 8.13 Let (uk) and (vk) be sequences in (IR ), s > 0, satisfying 
I151\} for some scale (fffc)- Suppose their Wigner measures at scale (fffc), /x 
and V are mutually singular. Then nrf' [uk + Vk] converges to + u. 

Proof. A proof of this result for s = may be found in or jl3j . 
For the general case, it suffices to take into account R,emark l8.4l to conclude 
that the Wigner measures of {£kDx)~^ Uk and {ekD^)"^ Vk are (O^^^A* 

^- These are clearly mutually singular and thus the aforementioned 
L^-version of the present result gives 

[{skD^r' {uk + Vk)] - (0"'^ ^ + (0"'' 



and finally 

m'" [uk + Vk] - (0'' ((0"'' + (0"'' u)=^Ji + v 
as claimed. ■ 

Acknowledgments: This article extends and improves some of the 
results of the author's Ph.D. Thesis. He wishes to thank the guidance of his 
Ph.D. advisor Enrique Zuazua. He also wishes to thank Patrick Gerard, for 
many helpful discussions and suggestions. 

This work has been supported by projects BFM02-03345 of MCyT (Spain) 
and HYKE (ref. HPRN-CT-2002-00282), HMS2000 (ref. HPRN-CT-2000- 
00109) of the European Union. 



52 



References 



[1] de Boor, C. A practical guide to splines. Revised edition. Applied Math- 
ematical Sciences, 27. Springer- Verlag, New York, 2001. 

[2] de Boor, C; DeVore, R.A.; Ron, A. Approximation from shift-invariant 
subspaces of (R°').. Trans. Amer. Math. Soc. 341(2), (1994), 787- 
806. 

[3] Brassart, M. The semi-classical limit in a crystal subject to exterior 
forces. Preprint (2003). 

[4] Evans, L.C. Weak convergence methods for nonlinear partial differential 
equations. CBMS Regional Conference Series in Mathematics, 74. Pub- 
lished for the Conference Board of the Mathematical Sciences, Wash- 
ington, DC; by the American Mathematical Society, Providence, RI, 
1990. 

[5] Gerard, P. Microlocal defect measures. Comm. Partial Differential 
Equations 16(11), (1991), 1761-1794. 

[6] Gerard, P. Mesures semi-classiques et ondes de Bloch. Seminaire sur 
les Equations aux Derivees Partielles, 1990-1991, Exp.No.XVI, Ecole 
Polytech., Palaiseau, 1991. 

[7] Gerard, P. Oscillations and concentration effects in scmilinear dispersive 
wave equations. J. Funct. Anal. 141(1), (1996), 60-98. 

[8] Gerard, P.; Leichtnam, E. Ergodic properties of eigenfunctions for the 
Dirichlct problem. Duke Math. J. 71(2), (1993), 559-607. 

[9] Gerard, P.; Markowich, P A.; Mauser, N.J.; Poupaud, F. Homogeniza- 
tion limits and Wigner transforms. Comm. Pure Appl. Math. 50(4), 
(1997), 323-379. 

[10] Hernandez, E.; Weiss, G. A first course on wavelets. Studies in Ad- 
vanced Mathematics. CRC Press, Boca Raton, FL, 1996. 

[11] Lions, P.-L. The concentration-compactness principle in the calculus 
of variations. The limit case. I. Rev. Mat. Iberoamericana 1(1), (1985), 
145-201. 

[12] Lions, P.-L. The concentration-compactness principle in the calculus of 
variations. The limit case. IL Rev. Mat. Iberoamericana 1(2), (1985), 
45-121. 



53 



[13] Lions, P.-L.; Paul, T. Sur les mesures de Wigner. Rev. Mat. Iberoamer- 
icana 9(3), (1993), 553-618. 

[14] Markowich, P. A.; Mauser, N. J.; Poupaud, F. A. Wigner-function ap- 
proach to (semi) classical limits: electrons in a periodic potential. J. 
Math. Phys. 35(3), (1994), 1066-1094. 

[15] Martinez, A. An introduction to semiclassical and microlocal analysis. 
Universitext. Springer- Verlag, New York, 2002. 

[16] Meyer, Y. Ondelettes et operateurs. I. Ondelettes. Actualites 
Mathematiques. Hermann, Paris, 1990. 

[17] Ron, A. Introduction to shift-invariant spaces. Linear independence. 
Multivariate approximation and applications, 112—151, Cambridge 
Univ. Press, Cambridge, 2001. 

[18] Tartar, L. if-measures, a new approach for studying homogenisation, 
oscillations and concentration effects in partial differential equations. 
Proc. Roy. Soc. Edinburgh Sect. A 115(3-4) (1990), 193-230. 

[19] Unser, M. Samphng-50 Years After Shannon, Proc. IEEE 88(4), 
(2000), 569-587. 

[20] Wigner, E.P. On the quantum correction for thermodynamic equilib- 
rium. Phys. Rev. 40, (1934), 742-759. 



54 



